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Abstract 

We discuss application of population synthesis for binary stars to pro- 
genitors of SNIa. We show that the only candidate systems able to sup- 
port the rate of SNe la ui a ~ 10 -3 yr _1 both in old and young populations 
are merging white dwarfs. In young populations (~ 1 Gyr) edge-lit deto- 
nations in semidetached systems with nondegenerate helium star donors 
are also able to support a similar ir a . The estimated current Galactic rate 
of SNIa with single-degenerate progenitors is ~ 10 -4 yr -1 . 



1 Introduction 

There is little doubt that explosions of SNIa are thermonuclear disruptions 
of the mass-accreting carbon-oxygen white dwarfs (CO WD) in binaries. The 
main facts arguing for this are: released energy per 1 g is comparable to eco— >Fej 
explosive nature of the events suggests that degeneracy plays a significant role; 
explosions may occur long after cessation of star formation; hydrogen is not 
detected in the spectra of SN la (but see discussion of single-degenerate scenario 
below) . 

Identification of the SN la progenitors is important for several reasons. It 
may help to constrain the theory of binary-star evolution. Modeling and un- 
derstanding the explosions requires knowledge of the initial conditions for them 
and the environments in which they take place. Evolution of the galaxies de- 
pends on the radiative, kinetic energy, and nucleosynthetic output of SNIa and 
the evolution of SN la rate in time, which, in turn, depend on the nature of the 
progenitor systems. The nature of the progenitors is related to the use of SNIa 
as distance indicators for determination of cosmological parameters Hq and go- 
Evolution of the luminosity function and the rate of SNe is important in this 
respect. 

A successful model for the population of progenitors of SN la has to explain 
the i nferred Galactic rate of events (4 ± 2) • 10~ 3 yr _1 l|CaPDellaro & Turattol 
1200 lh . the origin of the observ ational diversi ty among local (z < 0.1) SNela 
— 36 ± 9% may be "peculiar " ijhi et al.ll200llh and the occurrence of SNe la in 
stellar populations having a wide range of ages. 

Below, we discuss the scenarios of formation of binary systems in which SN la 
may occur and the rate of SNIa, predicted by different scenarios. 



2 Population synthesis 

The data provided by stellar evolution theory allows to construct numerical evo- 
lutionary scenario that describes the sequence of transformations of a binary sys- 
tem with given initial masses of components and their separation (M\q, M20, clq) 
that it can experience in its lifetime. 

Statistical studies of stars provide information on the binarity rate and the 
distributions of binaries over M10, ap, qo — M20/M10. Combined with star 
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Figure 1: Evolutionary scenarios for possible progenitors of SNIa. 



formation history, this allows to estimate the birthrate of the systems with 
a given set of (Mxo, M 2 o,ao) at any epoch. Then, it is possible to compute 
their contribution to the past or present population of stars of different types. 
Integration over whole space of initial parameters or Monte Carlo simulation 
for a large sample of initial "binaries" gives a complete model of the population 
of binaries and occurrence rates of different events, e. g., SN. Objects of the 
same type may be formed by several routes, hence, one may expect variations 
of SNIa. 
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Table 1: Occurrence rates of SNe la in candidate progenitor systems (in yr 1 ) 


Donor 


CO WD MS/SG He star He WD RG 


Counterpart 


Close Supersoft Blue AM CVn Symbiotic 
Binary WD XRS sd Star 


Mass transfer mode 


Merger RLOF RLOF RLOF Wind 


Young population 
Old population 

Young population 


Direct carbon ignition (Chandrasekhar SN) 
10 3 1CT 4 1CT 4 1CT 5 1(T 6 

10 3 icr 5 KT 6 

Edge-lit detonation (sub-Chandrasekhar SN) 

< icr 4 10 3 - < 10- 3 



3 Evolutionary scenarios for progenitors of SN la 

Figure ^ shows (not to scale) a simplified flowchart of the main scenarios in 
which one may expect formation of a progenitor of SNIa - a CO WD that 
may ignite carbon in the center. The rates of formation of potential SNIa via 
different channels are summarized in Tabled 

S cenario A ["double-degenerate"- DP - scenario, l|Tutukov fc Yungelsonl 

Il98li IWebbinklll98i llben fe Tutukovlll984l) ] starts with a main-sequence (MS) 
binary with Mio, M20 ~ (4—10) Mq. The system is wide enough for Roche-lobe 
overflow (RLOF) to occur when the primary is an AGB star with a degenerate 
CO core. After RLOF, a common envelope (CE) forms. If components do not 
merge inside CE, the core of the primary becomes a CO WD. After dispersal 
of CE, the system remains wide enough for the secondary to become a CO 
WD too. The angular momentum loss (AML) via gravitational waves radiation 
(GWR) results in the RLOF by the lighter of two WD. Mass loss proceeds on 
dynamical time scale and in several orbital revol utions Roche-lobe filling WD 
turns into a disk around the more massive WD ifTntukov fc YiingelsoTill97fll 



iBenz et 'all lTQQO'). If the total mass of the system exceeds Mch, accretion from 
the disk may result in accumulation of Mch by the "core" and SN la. 

Scenario B is realized in the systems with M20 ~ 2.5 Mq and such a sep- 
aration of components after formation of the first WD that the secondary fills 
its Roche lobe in the hydrogen-shell burning stage and becomes a helium WD. 
Like in scenario A, dwarfs are brought into contact by the AML via GWR. 
Unstable merger, most likely, results e ither in ignition of He at the interface of 
accretor and disk l|Ergma et aljEooi)) . forma tion of a CE and loss of He-ric h 
matter or in formation of an R CrB-type star l|Webbinklll98i llben et al Ell). 
If a stable semidetached system (of an AM CVn-type) forms, accumulation of 
Mch by accretor becomes possible. 

In scenario Scenario C ["edge-lit detonation" - ELD - scenario, l|Livnd 

Il990|) ] 2.5 M2O/M0 Zi 5 and the separation between components after the 
first CE phase is such that the secondary fills its Roche lobe before core He 
ignition and becomes a low-mass [~ (0.35 — 0.8) Mq] compact He-star. Low- 
mass helium remnants of stars have lifetime comparable to the MS-lifetimc of 
their progenitors. This allows AML via GWR to bring He-stars to RLOF be- 
fore exhaustion of He in the cores. If mass loss occurs stably, M a , ~ (2 — 3) ■ 
10~ 8 M<^ yr" 1 , almost independ ent of the mass of companion l|Savoniie et al.l 
Il986t iTutukov fc Fedoro^alll989|) . Under such M a a degenerate He-layer forms 
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Figure 2: Rates of potential SNIa-scale events after a 1-yr long star formation 
burst that produces 1 Mq of close binary stars. 



atop WD and detonates when its mass increases to ~ 0.1 Mq l|Limongi fc Tornambel 

Detonation of He produces an inward propagating pressure wave that 
leads to close-to-center detonation of C. The total mass of configuration in this 
case may be sub-Chandrasekhar. 

Scenario D ["single degenerate" - SD - scenario, llWhelan fc Ibedll97at> ] 
occurs in the systems where low-mass MS (or close to MS) stars [M20 ^ (2 — 
3) Mq] or (sub)giant (M20/M1 £ 0.8) companions to WD stably overflow Roche 
lobes. Accreted hydrogen burns into helium and then into CO-mixture. This 
allows to accumulate Afch- 

Scenario E is the only way to produce SNIa in a wide system, via accu- 
mulation of a He layer for ELD or Mc,h by acc retion of stellar wind matter in a 
symbiotic binary (|Tutukov fc Yungelsonlll97fih . 

Scenarios A - E are associated with binaries of different types and with 
different masses of components. This sets an "evolutionary clock" - the time 
delay between formation of a binary and SN la. Figure [21 shows the differential 
rates of SN la produced via channels A, C, and D after a burst of star formation. 
The DD-scenario is the only one that may operate in the populations of any 
age, while SD- or ELD-scenarios are not effective if star formation ceased several 
Gyr ago. 

Table 0] presents the order of magnitude model estimates for after 10 Gyr 
since beginning of star formation in the populations that have similar total mass 
comparable to the mass of the Galactic disk. Computations were made by the 
code used, e. g., by Tutukov and Yungelson (1994) and Yungelson and Livio 
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(1998) for the value of common envelope parameter a ce = 1. (Differences in 
assumptions in population synthesis codes or parameters of computations result 
in numbers that vary by a factor of several; this is the reason for giving only 
order of magnitude estimates). "Young" population had constant star formation 
rate for lOGyr; in the "old" one the same amount of gas was converted into 
stars in 1 Gyr. We also list in the table the types of observed systems associated 
with certain channel and the mode of mass transfer. Like Fig. [2 tabled shows 
that, say, for elliptical galaxies where star formation occurred in a burst, DD- 
scenario is the only one able to respond for occurrence of SN la, while in giant 
disk galaxies with continuing star formation another scenarios may contribute 
as well. 

For a certain time the apparent absence of observed DD with M to t > Mch 
merging in Hubble time was considered as the major "observational" difficulty 
for scenario A. Theoretical models predicted that it may be necessary to in- 
vestigate for bina rity up to 1000 field WD with V ,$ 16 -j- 17 for finding a 
proper candidate l|Nelemans et alJl200lh. The "nece ssary" number of WD was 
studied within SPY-proiect l|Napiwotzki et alJl2~00l|) and resulted in discovery 
of the first super-Chandrasekhar pair of dwarfs [R. Napiwotzki (this volume), 
Napiwotzki et al., (2003)]. 

On the "theoretical" side, it was shown for one-dimensional non-rotating 
models that the central C-ignition and SN la explosion are possible only for 
M a < (0.1 - 0.2)M Edd jNomoto fc IberJll985|) . But it was expected that in 
the merger products of bin ary dwarfs M a is close to MEdd ~ 10~ 5 M Q yr _1 
l|Mochkovitch & Liviol[i 990) because of high viscosity in the transition layer be- 
tween the core and the disk. For such M a the nuclear burning will start at the 
core edge, propagate inward and convert the dwar f into an ONeMg one. The 
latter will collapse without SNIa l|lsern et alJ ll983). However, consideration of 
the role of deposition of angular momentum into central object (Piersanti et 
al., 2003a, b) has shown that, as a result of spin- up of rotation of WD, insta- 
bilities associated with rotation, deformation of WD and angular momentum 
loss by distorted configuration via GWR, M a that is initially ~ 10 -5 M Q yr _1 , 
decreases to ~ 4 - 10 -7 M Q yr _1 . For this M a close-to-center ignition of carbon 
becomes possible. 

Because of long apparent absence of an observed "loaded gun" for the DD- 
scenario and its "theoretical problems", SD-scenario (D) is often considered 
as the most promising one. However, it also encounters severe problems. No 
hydrogen is observed in the spectra of SNIa, while it is expected that ^0.15 
Mm of H-rich matter may be stripped from the companion by the SN shell 
l|Marietta et al.112000^1 1 . Hydrogen may be discovered both i n very early and 
late optical spectra of SN and in rad io- and X-ray ranges llEck et all Il995l 
Mar ietta et alJl2000t iLentz et al.ll2002h . As well, no expected l)Marietta et alJ 
2000i ICanal et alJl2001r IPodsiadlow ski 2003) high luminosity and/or high ve- 
locity former companions to exploding WD were discovered as yet. 

In the SD-scenario, hydrogen first burns into helium and then into C/O mix- 
ture. However, two circumstances hamper accumulation of Afch- 

At M a 10~ 8 M Q yr _1 all accumulated mass is lost in Novae explosions l|Prialnik fc Kovet3 
Il995l) . Even if Af a allows accumulation of He-layer, most of the latter is lost 

1 Recently discovered SNIa 2001ic and similar 1997cy jHamuv^^al . 2003) may b elong to 
the so-called SN 1.5 type or occur in a symbiotic system JChugai h Yu ngclson 20041) . 
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Figure 3: The rate of accumulation of Met in the SD-scenario (in yr _1 ), de- 
pending on the masses of WD-accretors and MS- or SG-donors at the beginning 
of accretion stage. 

after He-flash l|lben fc Tutukovll 19961 fCassisi et al]ll99St IPiersanti et al.lll999(l . 

dynamically or via frictional interaction of binary components with giant-size 
CE. Thus, the results of computations strongly depend on the assumptions 
on the amount of mass loss in the nuclear-burning flashes. The flashes be- 
come less violent and more effective accumulation of m atter may occur if mass 
is transferred on the rate close to the therma l one l|lben fc Tutukovl 1 19841 
lYungelson fc Liviol Il998t llvanova fc Taaml l2003h . However, this assumption 
seems to lead to overprod uction of supersoft X- ray sources [see the estimate 
of the numb er of sources in iFed orova et alJ ll2004l) and completeness of surveys 
estimates in lDi Stefano fc Rappaportl l)l995|) ]. 

The "favorable" range of mass transfer rates w idens if mass exchan ge is sta- 
bilized by optically thick stellar wind from WD (Hac hisu et al.lll996|) . Under 
this assumption (not based on a rigorous treatment of the radiation transfer) , 
the excess of transferred matter over the upper limit for stable hydrogen burn- 
ing (~ 5 • 10 -7 M Q yr _1 for a 1M Q WD) is blown out of the system taking 
away specific angular momentum of the WD. This allows to avoid formation of 
CE for mass transfer rates up to ~ 10~ 4 M Q yr _1 and, simultaneously, implies 
stable hydrogen burning and reduces mass loss in helium burning flashes. Fig- 
ure shows the range of masses of donors and accretors in "successive" SN la 
progenitors at the beginning of accretion onto the WD stage, obtained under 
"stabilization" condition and for thermal-time scale mass transfer by Fedorova 
et al. (2004). The maximum of j/i a in the latter study is 2 ■ 10~ 4 yr _1 , i. e., it 
still does not exceed ~ 10% of the inferred Galactic v\ & . Han & Podsiadlowski 
(this volume) obtain for this channel the rate up to 1.1 • 10 _3 yr _1 , closer to the 
observational estimate. 

An important source for discrepant i\ & obtained for SD-scenario may be the 
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Figure 4: Comparison of efficiency of accumulation of matter by a 1 M white 
dwarf under different assumptions. See text for details. 

difference in the assumptions on the mass accumulation efficiency. As an ex- 
treme example, the upper panel of Fig. ^ shows the efficiency of accumulation 
of He and C+O if one takes into account stellar wind mass loss by dwarfs that 
burn hydrogen steadily, mass loss in Novae explosions after Prialnik and Kovetz 
(f995) and estimates of mass loss in helium flashes after Iben and Tutukov 
(f996); the lower panel shows efficiency of accumulation under prescriptions 
adapted by Han & Podsiadlowski [Fedorova et al. (2004) implemented an "in- 
termediate" case: assumptions on H-accumulation after Prialnik and Kovetz 
and assumptions on He-burning similar to Han & Podsiadlowski]. 

Scenario C may operate in populations where star formation have ceased 
no more than ~ 1 Gyr ago and produce SN at the rates that are comparable 
with the Galactic v\ a . But the outcome of ELD currently seems to be not 
compatible with observations of SNla. "By construction" of the model, the 
most rapidly moving products of explosions have to be He and Ni; this is not 
observed. The spectra produced by E LD are not compatib le with observations 
of the overwhelming majority of SNla llHoeflich et al]ll99ffl . On the theoretical 
side, it is possible that lifting effect of rotation that reduc es effective gravity and 
degeneracy in the helium layer may prevent detonation l)Laneer et al.ll2003|) . 

Channel B most probably gives a very minor contribution to the total SN la 
rate since typical total masses of the systems are well below Mch- 

The peculiarity of channel E is the behavior of M a from the wind: it is 
initially very low and grows, as companion to the WD expands. Typical initia l 
masses of WD in symbiotic stars are well below 1 M Q ijYuneelson et al.lll995(l . 
For them it is more likely to accumulate a helium layer that may be lost in a 
thermal flash than accumulate Mch- 
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4 Conclusion 



1. Only DD may secure the observed v\ a both in old and young populations. 
Merging pairs with Mi + M2 — Mch were discovered after search in a WD- 
sample of appropriate size. Account for effects of rotation may solve the problem 
of central ignition in the merger product. Crucially needed is a study of the 
physics of merger which follows development of shocks and turbulence in the 
"transition" zone, transfer of momentum, rotation effects upon evolution of the 
"core-disk" configuration. 

2. Edge-lit detonations in He-accreting systems can be responsible for SN la- 
scale events only in the populations younger than ~ 1 Gyr. Lifting effect of 
rotation may reduce the number and scale of ELD. 

3. Single-degenerate scenario may contribute a fraction (~ 10%) of all events 
in young or intermediate age populations. The major obstacle to SD-scenario 
are H and He thermal flashes. Predictions of the rate of SD-events have to be 
reconciled with the number of supersoft X-ray sources. A crucial test for SD- 
scenario would be detection of H which may be present due to the interaction 
of SN shell with companion or a "slow wind" of pre-SN. 

4. In the DD-scenario one may expect that exploding objects would differ in 
mass and central C-abundance. In SD-scenario all exploding WD most probably 
have Mch, but differ in central C. It is unclear whether these differences may 
explain the diversity of observed SN la. 

This study was supported by RFBR grant 03-02-16254 and Federal Program 
"Astronomy" . The author acknowledges financial support of IAU and JD5 SOC 
that enabled his participation in the meeting. 
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